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ABSTRACT 


This  paper  presents  an  overview  of  a  fracture  mechanics  approach  to 
some  of  the  most  frequently  encountered  matrix-dominated,  sub-laminate 
cracks  In  epoxy-based  composite  laminates.  By  f' sub-laminate'^  it  is  meant 
that  the  cracks  are  internal  to  the  laminate,  generally  invisible  macro- 
scopically;  but  are  much  larger  in  size  than  those  microcracks  considered 
in  the  realm  of  micromechanics.  The  origin  of  sub-laminate  cracks  is 
assumed  to  stem  from  the  coalesence  of  natural  material  flaws  (also  micro- 
cracks)  which  occur  under  a  certain  favorable  laminate  stress  condition. 
Thus,  the  modelling  of  the  mechanisms  of  sub-laminate  crack  initiation  and 
propagation  is  essentially  mechanistic  and  probabilistic  in  nature.  Some 
specific  results  from  several  analytical/experimental  investigations  using 
graphite-epoxy  laminates  are  presented  and  discussed  in  this  paper. 
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INTRODUCTION 


Failure  mechanisms  in  structural  composite  laminates  have  been  viewed 
at  several  dimensional  levels.  Consider,  for  instance,  the  curved  lamina¬ 
ted  panel  shown  in  Figure  1(a).  Failure  of  this  structured  component  may 
be  caused  by  a  loss  of  the  global  stiffness  when  the  applied  load  reaches 
a  certain  critical  value.  To  describe  analytically  the  associated  failure 
mechanisms  and  hence  to  determine  the  critical  load,  structural  mechanics 
methods  such  as  buckling  and  post-buckling  theories  are  employed,  which 
relate  the  change  of  the  structural  geometry  to  the  applied  load.  In  this 
type  of  analysis,  the  stiffness  of  the  laminated  panel  and  the  kinematics 
of  its  displacements  are  among  the  most  predominant  factors. 

On  the  other  hand,  the  same  panel  may  fall  due  to  a  loss  of  strength 
at  a  local  defect;  for  example,  at  a  bolt-hole.  Figure  1(b).  In  this  case, 
rupture  of  material  will  begin  at  the  hole  and  may  propagate  into  a  large 
crack  whenever  a  certain  load  condition  is  reached.  But  in  order  to  'define 
the  conditions  for  rupture  initiation  and  propagation,  a  knowledge  of  the 
actual  stress  field  around  the  hole,  and  the  physical  mechanisms  of  the 
material  rupture  process  is  essential.  To  this  end,  it  will  be  necessary 
to  focus  further  on  the  local  defect  area  at  a  much  smaller  dimensional 
level.  As  illustrated  in  Figure  1(c),  for  instance,  the  lamination  de¬ 
tails  of  the  panel,  such  as  layer  Interfaces  and  the  stacking  sequence 
are  now  identifiable.  Consequently,  quantities  of  size  comparable  to  the 
layer  thickness  become  important;  the  influence  of  an  interface  defect 
(i.e.  delamination),  a  translamlnar  crack,  a  fiber  split,  etc.  are  new  fac¬ 
tors  to  be  considered.  For  it  is  believed  that  the  development  of  rupture 
around  the  hole  is  precipitated  by  these  sub-laminate  cracking  activities. 


any  failure  analysis  performed  at  this  dimensional  level  muse  first  address 
the  individual  mechanisms  of  the  various  sub-laminate  cracks,  and,  then 
their  interaction  and  coalescence  mechanisms  when  they  occur  simultaneously 
and/or  sequentially. 

It  is  well  known  that  the  formation  and  propagation  mechanisms  of  a 
crack  are  governed  by  the  conditions  that  exist  in  the  close  vicinity  of 
the  crack  front.  In  the  case  of  the  sub-laminate  cracks  mentioned  above, 
further  focusing  of  the  crack  front  region  will  reveal  the  microstructure 
of  the  composite  system.  Figure  1(d).  Here,  a  clear  distinction  can  be 
made  between  the  fiber  and  the  matrix  phases  of  the  material.  At  this 
dimensional  scale,  one  also  finds  a  random  distribution  of  material  micro¬ 
flaws,  be  it  voids,  broken  fibers,  matrix  crack,  debonded  fiber-matrix  in¬ 
terfaces,  just  to  mention  a  few,  see  Figure  1(e). 

Although  the  inherent  microflaws  are  small  in  size,  usually  on  the 
order  of  the  fiber  diameter,  they  may  behave  as  local  stress  risers  under 
the  applied  load.  When  a  certain  local  condition  is  reached,  the  micro¬ 
flaws  Interact  and  coalesce  to  form  an  actual  crack,  much  larger  in  size, 
and  identifiable  at  the  sub-laminate  level.  Hence,  a  proper  analysis  of 
the  mechanisms  must  be  capable  of  delineating  the  individual  behavior  of 
the  various  types  of  stress  risers,  and  their  interactions  when  one  is  lo¬ 
cated  close  to  another.  All  these  depend  profoundly  on  the  probabilistic 
nature  of  the  microflaws  in  terms  of  their  size,  location  and  density 
distributions. 

The  physical  process  of  material  failure  as  portrayed  in  the  proceed¬ 
ing  discussions  is  seen  to  span  a  wide  range  of  dimensions.  Though  the 
entire  process  is  essentially  one  continuous  event,  failure  analysis  could 


be  conducted  only  within  Isolated  dimensional  regimes.  As  illustrated  in 
Figure  1,  the  entire  dimensional  spectrum  is  separated  roughly  into  three 
analysis  regimes;  namely  the  structural  mechanics ,  the  macromechanics  and 
the  micromechanics  regimes.  Within  each  analysis  regime,  the  material 
failure  process  Involves  a  distinct  set  of  material  and  geometrical  para¬ 
meters. 

Analytical  and  experimental  treatments  on  laminate  failure  at  the 
structural  mechanics  level  have  generally  been  approached,  using  the 
classical  lamination  theory  (see  [1]).  Stiffness-predominant  structural 
responses  such  as  panel  buckling  and  post-buckling  have  been  successfully 
analyzed  [2].  The  effect  of  local  material  damage  (e.g.  from  a  bolt-hole 
a  delamination  caused  by  impact,  etc.)  upon  the  global  response  can  also 
be  evaluated  [3].  But,  a  more  detailed  description  about  the  damage  form' 
ation  and  propagation  mechanisms  always  required  separate  treatment. 

For  a  class  of  epoxy-based  composite  laminates,  e.g.  graphite-epoxy 
systems,  material  damage  mechanisms  around  a  bolt-hole  or  near  the  tip  of 
a  line-notch  have  been  modeled  by  the  crack  growth  formula  that  are  de¬ 
rived  from  the  classical  fracture  mechanics  for  structural  metals  [4-6]. 
Essentially,  the  notched  laminate  is  regarded  as  a  2-dimensional  body, 
and  the  crack  growth  from  the  notch  is  assumed  self-similar  and  catastro¬ 
phic.  In  order  for  the  fracture  formula  to  correlate  with  a  body  of  ex¬ 
perimental  data,  the  notch  size  parameter  must  be  adjusted  to  include  an 
empirically  defined  "intense  energy  region"  at  the  crack  tip.  Thus,  to¬ 
gether  with  the  laminate’s  fracture  toughness  (Kc  or  Gc) ,  yet  to  be  deter 
mined  experimentally,  the  model  consists  of  two  disposable  parameters 
which  must  be  adjusted  to  fit  a  given  set  of  test  data. 


Clearly,  the  empiricism  in  this  approach  is  less  than  desirable.  As 
was  discussed  earlier,  the  development  of  material  damage,  say  around  a 
small  bolt-hole,  is  precipitated  by  a  multitude  of  sub-laminate  cracking 
activities  in  that  region.  The  term  "intense  energy  region"  is  in  fact  a 
gross  representation  of  these  activities  at  a  larger  dimensional  scale. 

It  would  be  ideal,  therefore,  to  carry  the  analysis  across  the  dimensional 
boundary  into  the  macromechanics  regime  (see  Figure  1) ,  so  as  to  under¬ 
stand  how  sub-laminate  crackings  lead  to  the  damage  development  around  the 
hole. 

Early  studies  on  fibrous  composites  dealt  mostly  failure  mechanisms 
in  the  micromechanics  regime.  For  Instance,  a  considerable  amount 
of  theoretical  and  experimental  treatments  was  given  on  the  subject  of  the 
fiber-matrix  interface  mechanics  [7  -  10].  Fracture  models  for  various 
microcracks  such  as  depicted  in  Figure  1(e)  are  characterized  in  terms  of 
energy  absorption  processes,  including  interface  debonding,  fiber  pull-out, 
matrix  crack  bridging  [11] ,  etc.  Although  the  modeling  of  these  various 
microflaws  is  basic  to  the  strengthening  mechanisms  of  the  composite  system 
at  the  fiber-matrix  level,  it  remains  to  Incorporate  these  "micro"  models 
into  some  analysis  which  includes  the  mutual  interaction  and  the  coales¬ 
cence  mechanisms  of  the  microflaws.  Indeed,  it  would  be  ideal,  again,  to 
have  a  fracture  model  which  accounts  for  both  the  mechanistic  failure 
behavior  of  the  microflaws  and  the  probabilistic  nature  of  their  existence; 
the  analysis  can  then  be  carried  into  the  regime  of  macromechanics. 

Several  mechanistic/probabilistic  failure  models  have  been  attempted 
in  the  past.  Zweben  and  Rosen  [12]  studied,  for  example,  the  tensile 
strength  characteristics  of  unidirectional  composites,  assuming  the  micro- 


flaws  a  distribution  of  local  broken  fibers.  The  basic  concept  can  be  ex¬ 
tended  to  multi-directional  laminated  systems,  but  the  complexity  in  geo¬ 
metry  as  well  as  the  multiplicity  of  failure  modes  make  this  type  of  effort 
so  far  unproductive. 

In  recent  years,  increased  attention  has  been  given  to  fracture  model¬ 
ing  within  the  general  confines  of  macromechanics,  as  illustrated  in  Fig¬ 
ure  1(c).  At  this  analysis  level,  the  individual  material  layers  in  the 
laminate  are  approximated  as  being  homogeneous  and  anisotropic,  where  no 
distinction  of  the  fiber  and  the  matrix  phases  need  be  made.  Essentially, 
the  stiffness  properties  of  the  material  layer  are  represented  by  some 
"effective"  values  characterized  experimentally  as  basic  material  constants, 
but  failure  in  the  material  layer  is  not  governed  by  a  set  of  strength  con¬ 
stants.  Rather,  it  is  determined  by  a  certain  fracture  event  which  occurs 
at  the  "sub-laminate"  level.  Since  by  a  homogeneous  approximation  the 
identity  of  the  Inherent  material  microflaws  also  becomes  lost,  only  their 
gross  effects  upon  the  formation  and  propagation  of  sub-laminate  cracks  can  be 
retained.  Thus,  an  empirically  defined  "effective"  crack  size  is  introduced. 

As  it  will  be  discussed  in  detail  later  in  this  paper,  the  macroscopic 
view  point  allows  a  rational  formulation  for  a  general  fracture  model  for 
a  class  of  sub-laminate  cracking  in  epoxy-based  composite  laminates. 

Early  studies,  notably  by  Corten  [13],  Wu  [14],  Kanninen,  et.  al.  [15], 
have  articulated  the  viability  of  such  an  approach.  Major  advances  have 
since  been  made*  due  mainly  to  the  ever-expanding  computational  capabilities 
and  ever- revealing  NDT  methods.  These  modern  facilities  have  provided  a 
means  for  a  mere  rational  correlation  between  experiments  and  analysis. 


In  a  series  of  papers  [16  -  22],  Wang,  Crossman,  et.  al.  developed  a 
unified  energy  model  within  the  context  of  macromechanics.  It  describes  the 


growth  mechanisms  in  a  class  of  matrix-predominant  sub-laminate  cracks. 

The  specific  cracks  considered  included  interlaminar  (delamination)  and 
translaminar  (fiber-split)  crackings,  found  most  prevalent  in  graphite- 
epoxy  laminates. 

This  paper  is  a  summary  of  their  major  findings  obtainec  ring  the 
course  of  several  analytical  and  experimental  investigations 

FRACTURE  CONSIDERATIONS  IN  MACROMECHANICS 
The  energy  model  developed  by  Wang  and  Crossman  is  formulated  on  the 
energy  release  rate  concept  of  the  linear  elastic  fracture  mechanics 
(LEFM) .  All  variables  in  the  model  are  defined  within  the  context  of  the 
macromechanics.  Since  fibrous  composites  possess  some  unique  characteris¬ 
tics,  special  considerations  had  to  be  given  to  defining  some  of  the  varia¬ 
bles. 

Application  of  the  Griffith  Criterion 

The  essence  of  the  LEFM  is  that  material  failure  is  not  defined  from 
the  stand-point  of  strength  as  a  constant  material  property.  Rather,  it 
is  determined  by  an  analysis  of  the  kinetics  of  the  actual  process  of 
fracture  propagation.  The  classical  result  of  Griffith  [23]  pertained  to 
an  elastic  plate  which  is  uniformly  stretched  in  one  direction  by  o.  The 
plate  has  a  through  crack  of  size  2a,  orientated  normal  to  the  direction 
of  o.  The  length  of  the  crack  is  assumed  small  but  finite;  and,  the 
material  is  homogeneous  and  isotropic.  Griffith  postulated  that  at  the 
instance  of  crack  extension,  a  loss  of  the  stored  elastic  strain  energy 
near  the  crack-tip  region  is  resulted;  this  energy  is  converted  into  the 


surface  energy  of  the  crack.  A  balance  of  energy  during  a  virtual  crack 
extension  leads  to  the  criterion 


3U  _  3S 
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(1) 


Equation  (1)  defines  the  general  condition  under  which  the  existing 
crack  begins  to  propagate  in  self-similar  mode. 

For  the  plate  problem,  the  critical  stress  at  the  instance  of  crack 
extension  can  be  derived  from  (1) ,  giving 


a 
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where  y  is  the  free  surface  energy  density  of  the  material. 

When  applied  to  engineering  problems,  the  Griffith  theory  is  often 
modified  for  practical  considerations.  For  instance.  Equation  (2)  becomes 
unbounded  as  a  ->  0.  Certainly,  no  real  material  can  sustain  an  infinite 

stress.  This  limitation,  however,  can  be  circumvented  by  invoking  the 

existence  of  material  flaws.  That  is  to  assume  for  the  material  some 

characteristic  distribution  of  flaws;  the  worst  of  which,  having  a  size 

of  2a  >  acts  like  a  real  crack.  It  then  determines  a  finite  critical 
o 

stress  according  to  Equation  (2) . 

Of  course,  flaws  do  exist  in  real  materials,  especially  in  fibrous 
composite  systems.  But  the  physical  identity  is  lost  at  the  dimensional 
level  where  the  analysis  is  performed.  The  quantity  aQ  can  be  defined 
only  empirically  as  an  intrinsic  material  property.  As  it  will  be  dis¬ 
cussed  later,  the  value  of  aQ  can  be  orders  of  magnitude  larger  than  the 
fiber  diameter  in,  say,  graphite-epoxy  systems. 


Another  practical  consideration  Is  related  to  the  definition  of  y, 
the  free  surface  energy  density  of  the  material.  For  crack  in  brittle 
material  such  as  glass,  then  y  is  as  defined.  For  most  other  engineering 
materials,  crack  extension  is  found  to  associate  a  certain  degree  of  in¬ 
elastic  deformation  near  the  crack-tip  region.  Furthermore,  the  crack 
extension  path,  or  the  crack  surface,  shows  a  certain  degree  of  rugged¬ 
ness,  depending  on  the  heterogeneity  of  the  material  viewed  at  a  micro¬ 
scale.  Early  studies  by  Irwin  [24],  Orowan  [25]  and  others  on  structural 
metals  considered  the  right-hand  side  of  Equation  (1)  the  irreversible 
work  required  to  create  a  unit  crack  surface  area.  Hence,  the  quantity  y 
can  be  interpreted  as  the  energy  dissipated  in  the  crack-tip  region  during 
crack  extension.  Clearly,  y  will  then  depend  on  the  inelasticity  as  well 
as  the  microscopic  heterogeneity  of  the  material  locally  (near  the  crack- 
tip)  .  It,  therefore,  must  also  be  regarded  as  an  intrinsic  property  of 
the  material  to  be  defined  at  the  macromechanics  level.  Conceivably,  y 
has  to  be  an  averaged  value  over  a  relatively  large  crack  surface  area 
for  fibrous  composites  in  order  for  it  to  be  a  material  constant. 

These  considerations  are  of  fundamental  importance  when  a  crack-like 
failure  is  modeled  at  the  macroscopic  dimensional  level.  For  only  in  this 
context  can  a  macroscopic  fracture  model  be  developed  along  the  rational 
arguments  of  the  classical  fracture  mechanics. 

In  common  practice,  the  right-hand  side  of  Equation  (1)  is  replaced 
by  Gc  (■  2y) ,  known  as  the  critical  energy  release  rate  of  the  material. 
The  left-hand  side  is  a  function  of  the  applied  load,  the  geometry  of  the 
body  and  the  size  of  the  crack.  Thus,  for  a  crack  undergoing  self-similar 
extension,  the  Griffith  criterion  is  expressed  as 


Accordingly,  Che  development  of  the  fracture  model  for  sub-laminate 

cracks  rests  upon  the  calculation  of  G(o,a)  and  the  physical  measurement 

of  G  .  The  latter  is  simply  a  material  constant, 
c 

The  Calculation  of  G(o,a) 

In  the  theory  of  the  LEEM,  the  singular  stress  field  near  a  crack  tip 
in  a  homogeneous,  isotropic  elastic  body  is  represented  by  analytical 
functions  in  the  theory  of  complex  variables  [26].  The  near-field  stresses 
are  obtained  for  three  particular  modes  of  the  crack  opening.  These  are 
known  as  the  opening  mode  (I) ,  the  slidding  mode  (II)  and  the  anti-plane 
shearing  mode  (III).  For  each  mode,  the  stresses  are  expressed  in  terms 
of  the  associated  stress  intensity  factor  K  [27];  and  consequently,  the 
associated  G(o,a)  is  computed  in  terms  of  K.  Since  the  relation  between 
K  and  G  is  one-to-one.  Equation  (3)  reduces  to  the  form  ^ 


K(a,a)  -  Kc  (4) 

Similar  relations  between  K  and  G  for  orthotropic  media  having  a  crack 
orientated  along  one  of  the  major  axes  can  also  be  obtained  [28].  But  the 
analytical  solutions  for  the  singular  stress  field  often  require  tedious 
mathematical  derivations. 

Direct  solution  methods  for  G(o,a)  have  been  available;  among  them 
are  the  well-known  J-integral  method  [29,  30]  and  the  method  of  virtual 
crack  closure  technique  by  Irwin  [31]. 

Irwin  [31]  observed  that  the  elastic  strain  energy  released  during  a 
virtual  crack  extension  Aa  is  equal  to  the  work  done  in  closing  it  again. 
The  inverse  problem  provides  the  solution  for  the  surface  tractions  a 


over  Aa.  The  crack-tip  energy  release  rate  is  then  represented  by 

S  ■  24lJo  (5'4“)da  (5) 

where  Au  is  the  crack  opening  displacement  vector  over  Aa. 

If  the  crack  extension  involves  all  three  modes  (I,  II,  III),  the 
vector  product  in  Equation  (5)  will  give  a  sum  of  three  scalars,  associated 
respectively  with  G^,  G  ,  and  G^j. 

The  virtual  crack-closure  representation  is  particularly  adaptive  to 
machine  computations.  Rybicki  and  Kanninen  [32]  suggested  a  2-dimensional 
finite  element  technique  to  evaluate  G  for  a  line  crack  in  a  plane.  The 
crack-tip  stress  vector  a  and  the  displacement  vector  Au  in  Equation  (5) 
are  approximated  by  the  nodal  forces  and  displacements  respectively,  in  a 
finite  element  representation  (for  detail,  see  [18]). 

Wang  and  Crossman  applied  this  technique  in  a  generalized  plane  strain 
finite  element  routine  [33],  which  can  s initiate  a  line  crack  propagation 
in  the  2-dimensional  cross-section  of  a  laminate.  Since  a  general  lami¬ 
nate  under  load  may  suffer  cross-sectional  warping,  the  routine  actually 
computes  3-dimensional  stresses  and  displacements  [33]. 

If  a  laminate  is  subjected  to  the  far-field  stress  oq  and  a  sub¬ 
laminate  crack  is  to  be  simulated,  the  crack-tip  energy  release  rate  can 
be  expressed  in  the  general  form 

G  ■  C  (a/t) *t* (a  /E  )2  (6) 

e  e  oo 

where  E  is  laminate  stiffness  in  a  direction,  and  t  is  the  linear  scale 
o  o 

between  the  actual  model  and  the  finite  element  model.  C  is  a  function 

e 

dependent  only  on  the  crack  size  a,  which  is  routinely  generated  for  a 


given  type  of  crack  in  a  given  laminate. 

Similarly,  if  the  laminate  is  subjected  to  a  uniform  temperature  load 
AT,  and  if  a  thermally  induced  crack  opening  is  resulted,  then  the  associa 
ted  energy  release  rate  at  the  crack  tip  is  expressed  by 

Gt  -  CT(a/t)-f(AT)2  (7) 

Generally,  the  laminate  is  prestressed  by  a  -AT  due  to  curing,  so  a 
combined  effect  is  resulted  when  <jq  is  applied; 

g  -  tCce)1/2*e0  +  (cT)1/2*AT]2*t  (8) 

where  e  is  the  far-f ield  strain  (*  a  /E  ) . 
o  o  o 

It  is  seen  that  the  finite  element  technique  is  extremely  versatile, 
and  can  be  efficiently  executed  for  simulating  complicated  sub-laminate 
crackings  such  as  delamination. 

Nonetheless,  the  accuracy  of  the  technique  has  been  a  subject  of  con¬ 
cern  of  many  analysts,  because  of  its  approximate  nature  in  representing 
a  mathematically  singular  stress  field.  As  has  been  demonstrated  by 
Raju  and  Crews  [34],  Spilker  and  Chou  [35],  and  Wang  and  Choi  [36],  the 
finite  element  stress  solutions  could  lose  significant  accuracy  in  the 
small  vicinity  of  the  singular  point;  the  region  in  which  the  stress  be¬ 
comes  inaccurate  is  generally  much  smaller  than,  say,  a  fiber  diameter, 
due  to  the  very  nature  of  the  singularity  [36].  However,  it  can  also 
be  demonstrated  that  the  stresses  in  that  small  region  do  not  contribute 
significantly  to  the  crack  opening  energy  release  rate,  especially  for  a 
crack  size  much  larger  in  proportion. 


Figure  2  shows  a  close  comparison  between  the  finite  element  computed 
G  and  the  elasticity  solution  counter-part,  for  a  transverse  crack  located 
In  the  mid-layer  of  a  3-layer  laminate  (see  Insert  of  Figure  2).  The 
extreme  layers  are  designated  as  material  1  and  the  mid-layer  as  material 
2.  Both  had  to  be  assumed  elastic  isotropic  materials.  The  exact  solu¬ 
tion  for  G  as  a  function  of  the  crack  size  a  is  given  by  Isida  [37],  using 
complex  stress  potentials;  while  the  finite  element  solution  is  computed 
using  a  rather  coarse  constant-strain,  triangular  element  mesh.  It  is 
seen  that  the  finite  element  solution  for  G(a)  compares  well  with  the  exact 
solution.  To  obtain  accurate  stress  closest  to  the  crack  tip,  one  can  still 
resort  to  the  finite  element  method  using  either  a  finer  mesh  or  a  higher- 
order  element  formulation  [38].  But,  such  extreme  measures  are  often  un¬ 
necessary  for  the  computation  of  G. 

The  Evaluation  of  G 

c 

When  fracture  occurs  in  the  material,  the  energy  released  in  the  pro¬ 
cess  is  expected  to  depend  on  the  morphology  of  the  fracture  surface, 
which  must  be  examined  at  the  microscale.  Fibrous  composites  are  known  to 
possess  complex  fracture  surface  details,  even  in  matrix-predominant  cracks. 
The  observed  delamination  surfaces  in  graphite-epoxy  composites  show,  for 
instance,  a  considerable  raggedness  because  the  crack  must  pass  around  the 
reinforcing  fibers.  But,  at  the  macroscopic  level,  the  crack  surface 
details  are  "smoothed  out";  and  their  effects  are  reflected  in  the  measured 
quantity  of  the  fracture  resistance  G£. 

For  this  reason,  G  measured  for  some  matrix-predominant  fractures  in 

c 

epoxy-based  composites  has  been  found  to  depend  on  the  direction  of  frac¬ 
ture  propagation.  Cullen  [39]  and  Williams  [40]  considered  two  different 


cases  of  delamination  as  illustrated  in  Figure  3.  The  first  case  is  0°/0° 

delamination  in  which  the  crack  path  is  in  the  fiber  direction,  while  the 

second  is  90°/ 90°  delamination  where  the  crack  path  is  transverse  to  the 

fiber  direction.  The  experiments  were  performed  using  a  graphite-epoxy 

unidirectional  laminate,  subjected  to  mode-I  cracking  condition.  They 

found  that  the  microscopic  morphology  of  the  90°/90°  delamination  surface 

exhibited  considerably  more  raggedness  than  the  0°/0°  delamination  surface. 

This  resulted  in  marked  differences  for  the  measured  G  .  Note  that  these 

c 

two  fracture  events  occur  essentially  in  the  same  interface  when  viewed 
macroscopically.  Yet,  the  respective  values  can  differ  greatly  depend¬ 
ing  on  the  direction  of  the  crack  propagation. 

Hence,  when  the  Griffith  formula  (3)  is  applied  for  cracks  in  compo¬ 
sites,  the  term  G  requires  a  precise  qualification.  Similarly,  when  a 
c 

test  method  is  devised  to  measure  G  ,  it  is  also  necessary  to  consider  the 

c 

dimensional  and  directional  characteristics  of  the  measured  data. 

Mode-I  interlaminar  G^c>  A  commonly  used  test  method  to  determine 

interlaminar  GJc  is  the  splitting  cantiliver  beam.  Cullen  [39]  and  Wilkins 

[41]  have  used  this  method  to  determine  the  Interlaminar  G^  when  the  crack 

is  propagating  in  the  direction  of  the  fibers  (0°/0°  delamination).  For 

some  graphite-epoxy  systems,  they  found  that  G^c  at  room  temperature  is 

about  0.85  lb/in,  or  130  J/m  .  This  value  is  about  twice  the  G_  measured 

Ic 


for  pure  epoxy  resin.  Williams  [40]  used  a  compact  specimen  which  simulates 
roughly  a  90°/90°  delamination.  He  found,  for  the  same  material  system, 
a  GIc  value  of  1.3  lb/in,  about  three  times  that  of  the  pure  resin.  Williams 
explained  that  the  fracture  surface  in  his  specimens  showed  fiber  breakage 
as  well  as  fractured  epoxy  debris;  this  had  resulted  in  a  higher  value 


for  G^c  than  that  found  for  0°/0*  delamination  by  Cullen  [39]. 
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In  another  paper  by  Wilkins,  et.  al.  [42],  it  is  reported  that  G^c 

in  delamination  of  0°/90°  interface  is  also  higher  than  that  of  0°/0° 

delamination.  These  findings  reaffirm  the  directional  dependent  nature 

of  G„ 

Ic 

Mixed-mode  interlaminar  G^  u)c*  7516  ®Plittin8  cantiliver  beam 
method  has  also  been  used  in  mixed-mode  cracking  experiment.  In  this  case, 
it  is  necessary  to  apply  different  loads  at  the  top  and  the  bottom  parts 
of  the  split  so  as  to  create  both  an  opening  (I)  and  a  slidding  (II) 
action.  Vanderkley  [43]  and  Wilkins  [41]  conducted  tests  on  the  same 
graphite-epoxy  system  (used  for  their  mode-I  tests) ,  and  found  the  total 
energy  release  rate  G^  ■  G^  +  G^^  which  exhibited  a  strong  de¬ 

pendence  on  the  G^/G  ratio. 

This  phenomenon  is  not  uncommon  in  mixed-mode  fracture.  Similar 
observations  were  reported  for  brittle  metals  as  well  as  pure  epoxy  resins 
[44].  It  is,  perhaps,  more  pronounced  in  fibrous  composites.  Generally,  it 
is  thought  that  the  increased  fracture  resistance  is  the  result  of  excess¬ 
ive  matrix  yielding  under  shear,  as  well  as  crack-closure  friction  due  to 
the  slidding  action. 

In  an  experiment  on  double-notched  off-axis  unidirectional  graphite- 

epoxy  laminates,  Wang,  et.  al.  [45]  measured  the  mixed-mode  G^  n)c  38  a 

continuous  function  of  G__/GT,  see  Figure  4.  It  is  seen  that  G/T  is 

II  I  (I,  H)c 

monotonically  increasing  with  G^/G^. 

But  in  the  same  experiment,  it  shows  also  that  the  opening  part  G 

•LC 

remains  roughly  independent  of  G^/G^.  This  suggests  that  mixed-mode  crack 
is  essentially  controlled  by  mode-I.  Of  course,  such  a  suggestion  is  only 
temporary;  more  study  is  needed  to  fully  understand  the  true  nature  of  the 
mixed-mode  crack  mechanisms. 


MECHANICS  OF  TRANSVERSE 


Physical  Mechanisms  at  the  Macroscale 

Transverse  cracks  are  found  in  epoxy-based  laminates,  even  at  a  low 

loading  level.  When  viewed  at  the  macroscopic  scale,  the  cracking  action 

is  simply  a  sudden  separation  of  fibers  by  breaking  the  epoxy  bond.  To 

illustrate,  consider  as  an  example  a  [0/90]  type  laminate  such  as  shown 

s 

by  the  Insert  in  Figure  2.  There,  material  1  is  the  0°-layer  and  material 
2  is  the  90°-layer.  Under  the  far-field  tensile  load,  material  2  could 
suffer  multiple  transverse  cracks.  Generally,  the  sequence  of  events  is 
as  follows:  A  single  crack  forms  first  when  the  applied  load  reaches  a 
certain  critical  value,  which  defines  the  "onset"  of  the  events;  as  the 
applied  load  increases,  more  similar  cracks  are  formed.  If  there  is  no 
other  failure  mode  setting  in  at  high  load  (e.g.  0°-rupture,  delamination, 
etc.),  t£e  number  of  transverse  cracks  will  continue  increasing,  until  it 
reaches  a  saturation  density. 

Figure  5  shows  a  load-sequence  x-radlographs  taken  for  a  [0/90] 

s 

graphite-epoxy  laminate  under  ascending  tensile  load.  Transverse  cracks 
in  the  90°-layer  are  seen  to  form  in  increasing  numbers.  From  these  photo¬ 
graphs,  a  plot  of  crack  density  (e.g.  cracks  per  inch  of  specimen  length) 
versus  the  applied  load  can  be  obtained.  Figure  6  shows  a  family  of  such 

experimental  plots  for  a  series  of  [0/90  /0]  laminates,  n  -  1,  2,  3,  4. 

n 

Examination  of  the  plots  in  Figure  6  reveals  several  interesting  fea¬ 
tures.  First,  the  onset  load  for  the  first  transverse  crack  seems  to  be 
Influenced  profoundly  by  the  90° -layer  thickness,  especially  when  it  is 
very  thin.  Take,  for  example,  the  case  of  n  ■  1  in  Figure  6,  the  onset 
load  is  almost  twice  that  of  n  ■  2. 


Secondly,  the  crack  density  also  shows  dependence  on  the  same  thickness 
factor.  Generally,  the  laminate  with  thinner  90°-layer  is  capable  of  yield¬ 
ing  a  higher  crack  density.  But,  for  the  case  of  n  ■  1,  failure  of  0°-layer 
at  high  load  interrupted  the  development  of  more  transverse  cracks. 

The  90°-layer  thickness  effect  on  transverse  cracking  was  first  docu¬ 
mented  experimentally  by  Bader,  et.  al.  [46].  They  attributed  the  effect 
to  the  constraining  actions  of  the  adjacent  0° -layers.  Observing  that  a 
transverse  crack  can  be  no  larger  than  the  90°-layer  thickness,  the  energy 
release  rate  which  drives  the  crack  is  thus  limited  by  the  same  factor. 

It  is  the  total  strain  energy  trapped  in  the  90°-layer  which  determines  the 
onset  of  the  cracking,  not  the  in-situ  tensile  stress. 

As  for  the  effect  on  crack  density,  it  has  been  explained  by  the  exist¬ 
ence  of  a  "shear-lag"  zone  at  the  transverse  crack  root  [46  -  48] .  That  is 
an  interlaminar  shear  stress  is  developed  on  the  0/90  interface  where  a 
transverse  crack  terminates.  This  shear  stress  is  singular  at  the  crack 
root,  but  decays  exponentially  a  distance  away  [20].  Similarly,  the  ln-situ 
tensile  stress  in  the  90 "-layer  is  nil  at  the  transverse  crack,  but  it  re¬ 
gains  its  far-field  value  outside  the  shear-lag  zone.  Thus,  ideally,  any 
two  adjacent  cracks  should  be  spaced  by  the  shear-lag  distance.  Since 
this  distance  is  proportional  to  the  90°-layer  thickness  [20] ,  hence  the 
observed  thickness  effect  on  crack  density. 

Although  the  shear-lag  concept  is  ideally  correct,  the  so-called 
"characteristic"  spacing  of  transverse  cracks  does  not  occur  in  practice. 
Often,  transverse  cracking  leads  to  other  failure  modes,  and/or  vice 
versa  [21]. 


AC  Che  microscopic  level,  Che  mechanisms  of  a  transverse  crack  are 
much  more  perplexed.  For  example,  Che  exacc  kinematics  of  crack  formation 
is  not  fully  known.  Post-tesc  SEM  examination  of  the  crack  surfaces  gen¬ 
erally  gives  a  very  ragged  appearance;  tiny  epoxy  debris  and  sometimes 
broken  fibers  are  also  seen  [40].  Figure  7  shows  an  x-ray  plane  view  of  a 

[+25/90.]  laminate  after  transverse  cracking,  left,  and  an  edge-view 
—  2  s 

micrograph  on  the  right.  It  is  seen  that  the  transverse  crack  is  practi¬ 
cally  a  plane  crack  of  a  rectangular  dimension,  which  is  bounded  by  the 
width  of  the  specimen  and  the  90°-layer  thickness.  There  is  no  evidence 
to  indicate  that  the  crack  formation  is  progressive  in  nature.  In  fact, 
all  experiments  tend  to  suggest  a  sudden  dynamical  formation. 

This  dynamical  nature,  though  still  conjectured,  has  been  qualified 
by  many  who  monitored  the  acoustics  emitted  during  the  crack  formation 
(see,  e.g.  [49]). 

While  it  is  difficult  to  reduce  these  physically  observed  facts  into 
a  general  law,  they  nevertheless  provide  the  necessary  rational  for  the 
formation  of  an  analytical  model.  In  what  follows,  an  energy  formulation 
is  presented,  which  describes  the  most  essential  observed  characteristics 
of  the  transverse  cracking  phenomenon. 

The  Energy  Formulation 

For  purpose  of  clarity,  consider  a  [0/90]  type  laminate,  as  shown  In 
Figure  8.  It  will  be  assumed  that,  in  the  90°-layer,  the  material  has  a 
random  distribution  of  microflaws.  The  gross  effects  of  the  aicroflaws  at 
the  macroscopic  scale  are  represented  by  a  characteristic  distribution  of 
"effective"  flaws  which  cannot  be  physically  seen.  But,  under  stress. 


these  effective  flaws  are  capable  of  propagating  suddenly  into  transverse 


cracks,  which  are  physically  real.  The  individual  size  of  the  effective 
flaws  is  denoted  by  2a;  and  any  two  adjacent  flaws  are  spaced  by  a  dis¬ 
tance  S,  see  Figure  8. 

Given  a  unit  length  of  the  specimen,  there  is  a  probability  density 
function  f(a),  and  a  probability  density  function  f(S).  For  sake  of  no 
evidence  to  suggest  otherwise,  the  two  functions  are  assumed  to  take  a 
form  for  normal  distributions  [50]; 

£<a)  ■  ITS  “»  [-<a  *  V2'20*2!  (9) 

f(s)  ■  sTS  [-<s  -  “s)2/2os2!  (10) 

where  u  and  a  are  the  mean  and  the  standard  deviation  of  the  respective 
distribution  functions. 

Among  the  effective  flaws,  the  size  of  the  "worst"  one  is  denoted  by 

2a  .  For  definiteness,  assume  a  the  99  percentile  of  f(a).  That  is  99Z 
o  o 

of  the  flaws  are  smaller  than  2a  . 

o 

Then,  under  the  far-field  load,  say  .e  ,  the  "worst"  flaw  2a  will 

o  o 

become  the  first  transverse  crack.  The  critical  value  for  e  at  the  onset 

o 

is  calculated  by  substituting  Equation  (8)  into  the  Griffith  criterion  (3) 
[(Ce)1/2.  eQ  +  (Op)172.  AT]2*  t  -  Gc  (11) 

where  C  and  C_  are  evaluated  at  a  ■  a  . 
e  T  o 

As  has  been  detailed  earlier,  the  energy  release  rate  coefficients 
Ce  and  CT  are  generated  numerically  by  the  finite  element  crack-closure 
routine,  given  the  geometry  and  material  moduli  of  the  laminas.  Thus, 
in  Equation  (11),  all  quantities  except  e  will  be  given:  AT  is  the  tern- 


perature  load  due  Co  cooling;  t  Is  Che  linear  scale  between  che  finice 


element  model  and  Che  accual  model;  and  Cakes  Che  value  of  GIc  which  is 

measured  for  mode-I  90°/90°  delaminaCion. 

Now,  che  difficulcy  resCs  upon  Che  choice  of  aQ,  or  for  Che  same  matter, 

Che  choice  of  p  and  o  in  che  discribucion  funcCion  (9). 
a  a 

Clearly,  2a^  muse  be  smaller  chan  che  90°-layer.  Buc  Che  laccer  can 
be  made  arbicrarily  large.  Hence,  a  finice  bound  on  aQ  exlscs  even  if  Che 
90°-layer  thickness  is  unbounded. 

As  an  example,  consider  Che  experimental  results  reported  in  [21] . 

The  censile  strength  of  T300/934  [90o]  laminates  averaged  a  -  7000  psi. 

If  the  same  "effective"  flaw  concept  is  assumed,  Chen  che  "worst"  flaw 
in  che  laminate  (parallel  Co  fibers)  determines  Che  strength  according 
to  che  Griffith  formula  (2),  yielding 

a  ■  G  E/ to  2  (12) 

o  c  u 

with  E  ■  1.7  X  10^  psi  and  G  ■  0.9  lb/ in  for  0°/0°  mode-1  delaminaCion, 

c 

Equation  (12)  determines  ag  -  0.01",  or  about  2  times  Che  ply  thickness  of 
the  T300/934  systems. 

It  may  be  assumed  that  for  Che  90°/90°  mode-I  delaminaCion  in  an  un¬ 
bounded  904-layer,  Che  "worst"  effective  flaw  size  should  be  no  greater 
Chan  aQ  .  0.01".  Indeed,  in  several  experimental  correlations  conducted 
by  Wang  and  Crossman  [19,  20,  21],  using  che  same  material  system,  aQ  is 
found  in  che  order  of  0.0075",  or  1.5  times  Che  ply  thickness. 

Numerical  Examples  by  Monte-Carlo  Simulation 

In  what  follows,  the  transverse  cracking  phenomenon  will  be  simulated 
numerically  by  the  so-called  Monte-Carlo  random  search  technic  ,.  .Tie 


considered  laminates  are  in  the  form  [0-/90  ]  ,  n  *  1,  2  and  3.  The 

2  ns 

material  system  is  the  AS-3501-06  graphite-epoxy  system  (for  the  material 
moduli,  see  [51]).  The  nominal  ply  thickness  of  this  system  is  0.0052". 

In  order  to  define  the  parameters  y  and  a  in  the  f(a)  and  f(S)  distri 
bution  functions,  the  following  values  are  chosen  for  the  laminates  n  ■  2 
and  3: 

y  -  0.0036",  a  -  0.0013"  (13) 

a  a 

y  -  0.0125",  a  -  0.0046"  (14) 

s  s 

These  are  chosen  to  reflect  the  fact  that  a  ,  being  over  99%  of  f(a), 

o 

may  take  a  value  (see  [50]), 

a  -  y  +  3a  _  0.0075"  (15) 

o  a  a 

and  that  for  yg  ■  0.0125"  it  Implies  80  effective  flaws  to  the  inch.  The 
choice  of  the  standard  deviations  is  a  matter  of  adjusting  to  the  actual 
data  scatter  [51]. 

As  for  the  laminates  of  n  ■  1,  the  thickness  of  the  90°-layer  is  only 

0.0104",  or  a  <  0.0052".  For  this  case,  the  choice  of  y  and  a  are  as 
o  a  a 

follows: 

y  -  0.0021",  a  -  0.0075"  (16) 

a  a 


so  that 


■A 


3 


a 


y  +  3c  0.0044" 


(17) 


The  choice  of  Che  spacing  parameters  remains  Che  same  for  n  “  1. 

The  Monte-Carlo  simulation  procedures  start  with  the  generation  of  a 
set  of  N  (  ■  80)  flaws  whose  random  sizes  and  spacings  are  represented  by 
the  respective  density  functions  (9)  and  (10).  This  is  done  by  generating 
first  a  set  of  N  random  values  in  the  interval  (0,1).  Then,  by  equating 
the  cumulative  function  of  f(a)  to  each  of  the  random  values,  a  random 
set  {a^},  i  *  1,  N  is  computed.  Among  the  values  in  {a^},  the  probability 
of  the  largest  to  be  equal  or  greater  than  a^  as  defined  in  (15)  or  (17) 
is  about  99%. 

Similarly,  a  random  set  {S^},  i  ■  1,  (N  -  1) ,  is  also  generated  by 
the  random  number  scheme.  is  then  assigned  to  be  the  spacing  between 
the  ith  and  the  (i  +  1)  th  flaws. 

Thus,  a  computer  research  follows,  which  determines  the  flaw  most 
likely  to  become  a  transverse  crack.  The  first  to  occur,  clearly,  is  the 
worst  flaw  in  the  {a^}  set,  and  the  corresponding  applied  load  for  the 
onset  of  the  first  crack  is  then  determined  using  Equation  (11) . 

After  the  first  crack  is  formed,  another  flaw  in  the  {a^}  set  will  be 

come  a  crack  at  a  slightly  higher  load.  But  this  flaw  is  located  at  a 

random  distance  from  the  first  crack.  And  the  presence  of  the  physical 

(real)  crack  has  a  stress  reducing  effect  on  the  rest  of  the  flaws.  Thus, 

if  the  second  flaw  to  become  a  crack  is  located  outside  the  shear-lag 

zone,  it  will  not  feel  the  presence  of  the  first  crack;  and  hence,  its 

available  energy  release  rate  at  the  instant  of  cracking  is  given  by 

G  (a, a),  the  same  as  in  Equation  (8). 
o 

On  the  other  hand,  if  the  second  flaw  is  located  inside  the  shear-lag 
zone,  then  the  available  energy  release  rate  is  reduced  by  a  factor  de- 


pending  on  Its  distance  from  the  first  crack. 

Generally,  for  the  flow  to  become  the  kth  crack,  one  must  search  to 
the  left  and  to  the  right  for  the  nearest  cracks.  If  so,  the  energy  release 
rate  at  the  Instance  of  the  kth  crack  is  given  by 

Gk(a,a)  -  R(S£)  GQ(a,a)  R(Sr>  (18) 

where  S0  is  the  distance  to  the  left  crack,  S  the  distance  to  the  right 
crack. 

In  Equation  (18),  R(S)  is  the  energy  rate  retention  factor  which  is 
generated  by  the  finite  element  routine  for  a  flaw  of  2a  (<2a^)  and  is 
placed  at  varying  distance  S  from  a  transverse  crack,  see  Figure  9. 

For  the  laminate  family  [0-/90  ]  ,  a  single  R(S)  curve  can  be  genera- 
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ted  if  expressed  in  terms  of  S/nt,  t  being  the  thickness  of  one  90“-ply, 
see  Figure  9.  Note  that  100%  energy  rate  retention  is  expected  beyond 
S  ■  9nt;  the  latter  is  actually  the  size  of  the  shear-lag  zone  [19]. 

The  applied  load,  corresponding  to  the  kth  crack  is  determined  from 

Gfc(a,a)  -  Gc  (19) 

The  search  essentially  simulates  the  cracking  process  as  it  would 
occur  naturally.  Each  random  simulation  represents  an  actual  test  case*, 
and  repeated  simulations  present  actual  tests  on  replicas. 

Figures  10,  11,  and  12  show  the  simulated  crack-density  versus  load 
plots  for  respectively,  811(1  ^^O^]  In  the 

simulations,  the  quantities  appearing  in  Equation  (11)  are  assigned  as 


follows: 


1.3  lb/in 


(20) 


AT  -  225°F,  t  -  0.0052",  G  - 

c 

The  shaded  band  In  each  of  these  plots  is  the  corresponding  experi¬ 
mental  data  band  from  4  test  specimens.  These  test  data  were  reported 
earlier  in  Reference  [52];  and  the  details  for  the  simulation  computer 
routine  are  reported  in  Reference  [51]. 

MECHANICS  OF  FREE  EDGE  DELAMINATION 
The  Classical  Free-Edge  Problem 

The  free  edge  delamlnatlon  problem  has  attracted  increased  research 
interests  ever  since  the  advent  of  composite  laminates.  The  phenomena  are 
frequently  observed  as  the  most  damaging  sub-laminate  failure  mode.  Gen¬ 
erally,  it  is  a  plane  crack  which  forms  along  the  laminate  free  edge  and 
propagates  inward  along  an  interface  of  two  adjacent  layers. 

Figure  13  shows  an  x-ray  plane  view  of  a  [+45/0/90]  graphite-epoxy 
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laminate  under  uniaxial  tension  (left) .  Delamination  is  seen  to  occur 
along  both  edges  of  the  laminate,  with  essentially  uniform  growth  toward 
the  center.  An  edge  view  of  this  crack  is  shown  on  the  right,  which  indi¬ 
cated  that  the  crack  is  basically  contained  inside  the  90°-layer,  not 
necessarily  in  any  one  given  layer  interface.  Moreover,  the  cracked  plane 
is  quite  zig-zagged  along  the  length  of  the  free  edge. 

Similar  photographs,  taken  for  a  [+25/90.  ,_]  laminate  after  failure, 

—  1/ 1  s 

are  shown  in  Figure  14.  Here,  the  edge  crack  is  seen  to  have  formed  in¬ 
side  the  90°-layer,  but  branched  to  the  25/90  interface  as  it  propagated. 
The  branch-out  is  due  to  the  skewed  cracks  which  occur  in  the  90°-layer 
ahead  of  the  edge  delamination  front. 


Note  that  both  laminates  have  a  90°-layer;  and  they  could  suffer  trans 
verse  cracking  under  the  tensile  load.  But,  because  the  thickness  of  the 
90°-layer  is  thin,  a  high  load  is  required  for  transverse  cracking  accord¬ 
ing  to  the  energy  theory  presented  earlier.  Instead,  in  this  case,  edge 
delamination  is  induced  as  the  first  event  of  sub-laminate  failure. 

Analytical  studies  of  the  free  edge  delamination  problem  have  origina¬ 
ted  from  the  much  celebrated  work  of  Pipes  and  Pagano  [53] ,  who  formulated 
and  computed  the  boundary- layer  interlaminar  stress  solutions  for  a  long 
symmetric  laminate  under  tension.  The  formulation  is  based  on  the  macro¬ 
scopic  ply-elasticity  theory,  which  regards  the  material  layers  as  indi¬ 
vidually  homogeneous  media;  material  and  geometrical  discontinuities  exist 
only  across  the  layer  interfaces. 

Generally,  the  free  edge  stress  field  is  three-dimensional,  and  is 

singular  at  the  intersect  of  the  free-edge  with  the  layer  interface  [36]. 

Hence,  interface  delamination  is  caused  by  the  highly  concentrated  edge 

stresses,  especially  the  interlaminar  normal  stress  a  . 
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In  a  series  of  tensile  strength  tests,  Bjeletich,  et.  al.  [54] 

examined  the  failure  modes  of  six  families  of  quasi-isotropic  laminates  by 

alternating  the  stacking  sequence  of  the  0°,  90°  and  +45°  layers.  Edge 

stress  analysis  indicated  a  compressive  a  along  the  free  edge  of  the 
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[0/90/+45]  laminate,  while  a  tensile  a  for  the  [+45/0/90]  laminate.  The 
latter  developed  prematured  delamination;  and  the  growth  of  it  had  led  to 
a  much  lowered  tensile  strength.  Clearly,  knowledge  of  the  free-edge 
stresses  can  provide  an  explanation  why  delamination  occur;  but  a  quanti¬ 
tative  prediction  for  its  occurrence  requires  a  more  precisely  defined 


criterion. 


In  a  study  on  delamination  for  similar  graphite-epoxy  laminates, 
Rodini,  et.  al.  [55]  determined  experimentally  the  critical  tensile  loads 


at  the  onset  of  free  edge  delamination  in  a  [+45  /0/90  ]  ,  n  ■  1,  2,  3 

—  n  n  ns 

family.  They  found  that  the  critical  laminate  tensile  stress  varied 

greatly  with  the  value  of  n.  Specifically,  the  critical  decreases  at 

the  rate  of  about  /n,  even  though  an  edge  stress  analysis  yields  identical 

o  for  the  same  a  for  all  values  of  n. 
z  x 

Such  a  90°-layer  thickness  dependent  behavior  is  similar  to  that  found 
in  the  transverse  cracking  problems,  suggesting  the  observed  phenomenon 
is  again  fracture  in  nature. 

The  Energy  Criterion 

In  the  work  of  Wang  and  Crossman  [18],  it  is  assumed  that  material 
flaws  exist  randomly  on  any  one  of  the  interfaces  between  the  material 
layers.  Those  flaws  which  are  located  within  the  free-edge  stress  zone 
form  an  "effective"  flaw  having  a  size  aQ  at  the  instance  of  onset  of  de- 
lamination,  see  insert  in  Figure  15.  Following  the  finite  element  crack- 
closure  procedure  by  allowing  virtual  extension  of  aQ  along  the  layer  in¬ 
terface,  an  energy  release  rate  curve  is  generated,  such  as  shown  in 
Figure  15. 

It  is  seen  that  the  energy  release  rate  G  increases  sharply  with  a, 

but  reaches  an  asymptote  at  a  ■  a  .  Generally,  the  value  of  a  is  about 

m  m 

one-half  the  layer-thickness  which  contains  the  delamination  [18].  Thus, 

the  layer  thickness  affect  the  value  of  a  ,  and  also  the  value  of  G. 

m 

Expressing  the  computed  energy  release  rate.  Equation  (6), 

G(eQ,  a)  -  CgCa)  •  t  •  eQ2 


(21) 


and  applying  the  Griffith  criterion  (3)  for  the  onset  of  delamination,  one 
obtains  the  critical  far-field  strain. 

1/2 

(e  )  r  -  [G  /C  (a)  *  t]lf,  a  -  a  (22) 

o  cr  cl  o 

Note  that  in  (22),  the  thermal  residual  effect  is  not  included.  In 
most  laminates,  thermal  effect  on  delamination  is  minimal  [56]. 

The  problem  comes  down  to  two  important  questions;  namely,  what  value 
is  aQ,  and  in  which  interface  is  it  located? 

To  answer  the  first  question,  recall  the  earlier  discussions  about  aQ 
in  unbounded  90°  laminates.  As  inferred  from  the  Griffith's  criterion,  a 

o 

is  in  the  order  of  0.01",  or  2  times  the  ply  thickness  for  commercial 
graphite-epoxy  systems.  It  is  believed  that  aQ  for  delamination  is  at 
least  of  this  magnitude  if  not  larger,  because  of  possible  additional  cut¬ 
ting  flaws  along  the  free  edge.  In  any  event,  if  am  in  Figure  15  is  less 
than  3  times  of  the  ply  thickness,  one  simply  uses  a  -  a^  and  predicts 
from  Equation  (22)  the  minimum  possible  load  for  the  onset  of  delamina¬ 
tion.  Thus 

<eo>cr  ’  “VW  '  t]1/2  (23) 

The  answer  to  the  second  question,  however,  is  complicated,  if  not 
Inconclusive.  Consider  as  an  example  a  [90^/0^/+< 452]  under  uniform  compres¬ 
sion.  Interlaminar  tensile  a  is  developed  in  this  laminate  and  edge  de- 
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lamination  is  induced  before  global  buckling,  see  Figure  16.  A  3-dimen- 

sional  through-thickness  display  of  o  is  shown  in  Figure  17.  It  is  seen 

z 

that  the  largest  o z  is  located  on  the  45/-45  interface  (in  fact,  this  <?z 
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is  singular  at  the  free  edge) •  On  the  other  hand»  o  along  the  mid-plane 

(z  ■  0)  is  finite  at  the  free  edge.  In  addition,  there  is  also  a  singular 

interlaminar  shear  stress  t  along  the  45/-45  interface  as  the  most  probable 

xz 

analysis  suggests  the  0/45  interface  as  delamination  site. 

However,  according  to  the  energy  release  rate  argument  as  well  as 
the  "effective"  flaw  hypothesis,  one  must  reply  upon  the  calculated  energy 
release  rate  curve  such  as  shown  in  Figure  15,  in  order  to  define  which 
one  of  the  interfaces  is  likely  to  delaminate. 

The  finite  element  crack-closure  results  of  the  G(a)  curves  for  this 
example  problem  is  shown  in  Figure  17.  Indeed,  the  interface  that  yields 
the  highest  energy  release  rate  is  th«  mid-plane  (z  ■  0) ,  not  the  45/-45 
interface.  Besides,  a  delamination  in  the  mid-plane  is  in  mode-I,  while  it  is 
mixed-mode  (I,  III)  in  the  45/-45  interface.  As  has  been  discussed  earlier, 
mixed-mode  is  usually  higher  than  'G^ .  Thus ,  the  mid-plane  is  the  pre¬ 
dicted  delamination  site*  not  the  0/ 45  interface. 

This  conclusion  is  confirmed  by  experiments  conducted  by  Wang  and 
Slomiana  [57],  who  also  investigated  compression-induced  delamination  in 
[09/90-/+45-]  and  [0/90/0/ 90/+45/+45]  laminates.  In  all  cases,  mid- 
plane  and  mode-I  delamination  was  predicted  and  observed.  Using  the  energy 

criterion  of  Equation  (23),  and  setting  G  ■  1.3  lb/in,  their  predicted 
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minimum  onset  loads  (o^)  compared  well  with  the  experimental  findings 
(average  of  3  specimens) : 


Laminates 


[902/02/±452]3 
[02/902/+452]g 
[0/ 90/ 0/90/+45/+45 ] g 


Predicted  onset  of  o 


43.9  ksi 


Experimental  Findin 


45.7  ksi 


In  the  previous  examples,  edge  delamination  is  the  first  and  only  sub- 
laminate  failure  mode  before  laminate  buckling.  The  onset  load  for  delami¬ 


nation  is  synonymous  to  the  final  laminate  failure  load.  In  most  other 
cases,  edge  delamination  and  transverse  cracking  are  often  two  competing 
failure  modes,  which  interact  each  other  through  a  complicated  mechanisms. 

Examine  again  the  photomicrographs  shown  in  Figures  13  and  14,  where 
a  [+45/0/90 ]g  and  a  [+25/90^2^3  laminates  were  loaded  by  axial  tension. 

The  thickness  of  the  90°-layer  in  both  laminates  is  so  small  that  trans¬ 
verse  cracking  is  possible  only  at  a  high  load.  Instead,  an  edge  delami¬ 
nation  is  induced  as  the  first  sub-laminate  failure  mode.  Nonetheless,  the 
stress  field  in  the  90°-layer,  which  contains  the  delamination,  appears  to 
be  extremely  complex.  The  skewed  cracks  shown  in  Figure  14  are  indicative 
of  the  complex  stress  field  ahead  of  the  delamination  front.  But,  just 
how  profoundly  these  secondary  cracking  events  influence  the  delamination 
growth  is  far  from  known. 

On  the  other  hand.  Figure  18a  shows  an  x-ray  plane  view  of  a  [HK25/ 90^ ] 
laminate  after  delamination.  Actually,  this  laminate  suffered  first  trans¬ 
verse  cracks  in  the  90°-layer;  subsequent  edge  delamination  in  the  25/90 
interface  developed  at  a  higher  load.  It  is  seen  from  the  picture  that 
delamination  actually  caused  many  more  transverse  cracks  because  of  stress 
concentration  along  the  curved  delamination  front.  In  the  area  where  de¬ 
lamination  had  not  occurred,  only  a  few  transverse  cracks  are  visible. 

This  phenomenon  has  repeatedly  been  observed  in  experiments  [21];  but  no 
serious  account  has  been  attempted  to  analytically  modelling  it. 


Law  [56]  assumed  chat  when  transverse  cracks  occur  in  the  90°-layer 


the  tensile  modulus  E„  and  the  Poisson  ratio  v_T  of  the  908-layer  both 

T  TL 

reduce  to  nearly  zero.  With  reduced  moduli  in  the  90°-layer,  the  overall 
laminate  becomes  energetically  delamination  prone  in  the  25/90  interface, 
resulting  in  a  mixed-mode  cracking  (I,  II).  Law  was  able  to  predict  the 


onset  loads  for  this  and  other  similar  cases  in  an  approximate  manner. 

Figure  18b  shows  an  x-ray  plane  view  for  a  [+25/90,]  laminate  soon 
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after  transverse  cracking  took  place.  In  this  case,  scattered  transverse 

cracks  appeared  first,  which  immediately  caused  delamination  on  both  25/90 

interfaces  in  the  area  where  they  occurred.  The  delamination  is,  in  fact, 

a  combined  effect  stemming  from  the  transverse  crack-tip  and  free  edge 

stress  concentrations.  Figure  19  illustrates  this  combined  effect  in  an 

isometric  view.  The  interface  shear  stress  t  ,  which  is  developed  at  the 
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transverse  crack  root,  tends  to  delaminate  the  25/90  interface  in  the  x- 

direction;  and  the  interlaminar  normal  stress,  a  ,  which  is  a  free  edge 
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effect,  tends  to  open  up  the  same  interface  in  the  y-direction.  The  com¬ 
bined  action  gives  rise  to  a  mixed-mode,  delta-shaped  delamination.  Evi¬ 
dence  of  this  type  of  delamination  can  be  found  by  examining  Figure  18b. 

Law  [56]  approximated  the  combined  effect  by  first  computing  for 
the  x-directlon  delamination  by  a  2-dimensional  analysis,  and  then  by 
computing  for  the  y-direction  as  in  the  case  of  free  edge  delamination. 

Since  both  G  and  G  are  considered  crack-driving  forces,  their  vector 
x  y 

sum  is  used  to  predict  the  onset  of  the  delta-shaped  delamination.  Ob¬ 
viously,  such  an  approximate  approach  is  only  tentative;  the  problem 
remains  to  be  treated  more  rigorously  by  a  3-dimensional  analysis  which 
can  simulate  a  contoured  plane  delamination.  But,  as  it  will  be  shown  in 


the  next  section.  Law's  approach  was  at  least  qualitatively  correct. 


It  is  now  apparent  that  sub- laminate  crackings,  especially  delamination 


are  physically  complicated  even  under  simple  tension.  The  energy  model 
presented  earlier  may  be  applied  only  in  cases  where  a  single  failure  mode 
is  present.  Interactions  amongst  multiple  failure  modes  are  yet  to  be 
investigated  more  extensively. 

An  Experimental  Case  Study 

An  experimental  case  study  on  transverse  cracking,  free  edge  delamina¬ 
tion  and  their  mutual  interaction  was  conducted  by  Crossman  and  Wang  [21]. 
They  chose  a  family  of  graphite-epoxy  laminates  (T300/934)  in  the  form  of 

[+25/90  ]  .  By  varying  the  single  parameter  n  -  1/2,  1,  2,  3,  4,  6  and  8, 
—  ns 

they  were  able  to  control  the  occurrence  of  the  various  competing  failure 
modes  in  the  laminates  when  loaded  in  tension.  Both  x-radiography  and 
microphotography  were  used  to  monitor  the  sub-laminate  cracking  development 
as  a  function  of  the  applied  load.  For  clarity,  a  graphical  summary  of 
their  experimental  findings  is  shown  in  Figure  20.  For  n  -  1/2  and  1,  the 
laminates  are  classified  as  having  a  ’’thin"  90°-layer.  Transverse  crack¬ 
ing  is  suppressed  due  to  the  lack  of  sufficient  strain  energy  stored  in 
the  90°-layer.  Hence,  free  edge  delamination  occurs  as  the  first  sub- 
laminate  failure  mode.  For  n  ■  1/2,  a  stable  mode-I  delamination  is  pro¬ 
duced  during  the  loading  until  final  failure.  For  n  *  1,  however,  trans¬ 
verse  cracks  are  Induced  in  the  area  of  delamination.  Because  of  these 
cracks,  the  delamination  plane  is  often  branched  to  the  25/90  interface. 
This  latter  appearance  is  indicative  of  an  interaction  taking  place  between 
the  two  failure  modes. 

When  n  ■  2  and  3,  the  laminates  are  denoted  as  having  "thick"  90°- 
layer.  In  both,  transverse  cracks  occur  first  long  before  delamination. 

It  may  be  thought  that  the  onset  of  transverse  cracking  is  an  independent 


event,  while  delamination  is  caused,  at  least  partially,  by  the  presence 
of  the  multiple  transverse  cracks.  Indeed,  the  delamination  plane  in  the 
case  of  n  *  2  resembles  that  for  n  *  1;  and  the  delamination  plane  in  the 
n  ■  3  laminate  is  entirely  in  the  25/90  interface.  Also,  for  n  •  3,  the 
growth  of  delamination  becomes  more  rapid. 

The  laminates  with  "very  thick"  908-layer  are  those  for  n  ■  4,  6, 
and  8.  In  all  cases  transverse  cracks  occur  first  followed  immediately 
by  delta-shaped  delaminations  emanating  from  the  transverse  crack  roots. 
These  soon  coalesce  to  form  large  scale  25/90  delamination,  resulting  in 
prematured  laminate  failure. 

Figure  21a  sunmarlzes  the  various  onset  loads  for  each  types  of  the 

laminates  (signified  by  the  value  of  n) .  Since  the  laminate  stress-strain 

response  for  all  cases  is  essentially  linear  until  final  failure,  it  is 

convenient  to  use  the  laminate  tensile  strain,  e  as  a  measure  of  the 
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applied  load.  Note  that  except  for  n  -  1/2  and  1,  transverse  cracking  is 

the  first  sub-laminate  failure  event.  Generally,  growth  of  multiple  cracks 

is  fully  developed,  followed  by  interlayer  delamination.  However,  as 

n  >  3,  the  load-gap  between  the  transverse  cracking  and  delamination  events 

grows  closer.  In  fact  for  n  ■  8,  these  two  events  are  so  closely  linked 

as  to  cause  Immediately  final  failure. 

Note  also  that  except  for  n  ■  1/2,  1  and  2,  onset  of  delamination  and 

final  failure  are  separated  only  by  less  than  5 X  of  loading. 

By  means  of  the  energy  model  presented  earlier.  Law  [56]  conducted  a 

computer  simulation  of  the  transverse  cracking  and  delamination  processes 

in  the  [+25/90  ]  family.  The  numerical  results  are  shown  in  Figure  21b, 

*  n  s 

which  is  superscribed  on  the  experimental  results  shown  in  Figure  21a. 
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Lav  has  conducted  four  different  types  of  crack  simulations:  (a)  onset 
of  transverse  cracking  in  the  90 “-layers  for  all  values  of  n;  (b)  mid-plane 
(z  ■  0) ,  mode-I  delamination  for  n  ■  1/2,  1,  2,  and  3;  (c)  25/90  interface 
delamination  for  n  ■  2,  3,  and  4  (the  90°-layer  is  saturated  with  trans¬ 
verse  cracks;  and  is  represented  by  reduced  moduli);  and  (d)  transverse 
crack  root/ free  edge  combined  delamination  of  the  25/90  interface  for  n  ■  3, 
4,  6  and  8  (as  illustrated  in  Figure  19). 

It  is  seen  from  Figure  21b  that  the  predicted  onset  loads  for  transverse 
cracking  agree  well  with  the  experiment,  while  the  predicted  delamination 
loads  (represented  by  3  curves)  form  an  envelope  to  which  the  experimental 
points  fall  below.  Although  the  predicted  delamination  loads  for  n  ■  6  and 
8  are  off  considerably  from  the  experimental  points,  the  energy  analysis 
seems  to  capture  the  physical  mechanisms  of  these  various  failure  modes; 
the  analysis  gives  correctly  the  trend  of  the  failure  loads. 

Final  Failure  Mechanisms 

In  the  case  study  discussed  above,  the  sub-laminate  failure  sequence 
is  seen  to  lead  to  the  final  rupture  of  the  laminate.  It  is  generally 
thought  that  final  fialure  is  determined  practically  by  the  strength  of 
the  load-carrying  plies;  i.e.  in  the  case  of  the  [+25/90  ]  family,  the 
+25-plies  should  carry  essentially  all  the  load.  But,  what  determines  the 
strength  of  the  +25-plies?  And,  is  it  adversely  influenced  by  the  sub¬ 
laminate  cracking  events  that  proceeded  the  final  laminate  rupture?  These 
questions  relate  to  the  final  failure  modes  and  their  associated  fracture 
mechanisms . 

Crossman  [58]  conducted  recently  an  extensive  experimental  investiga¬ 
tion  into  the  failure  sequence  of  several  families  of  laminates  under  uni- 


axial  tension.  Some  of  his  results  may  help  to  answer  the  above  questions. 

Consider  first  the  effect  of  thickness  of  the  load-carrying  plies  by 
examining  the  results  in  Table  A. 

It  is  noted  that  the  final  failure  mode  of  [+25/+25]g  is  due  simply 

to  fiber  breaking  across  the  laminate,  while  the  final  failure  mode  of 

[2S-/-25-]  involves  a  major  transverse  crack  in  the  center  core  (-25,°- 
2  2s  4 

layer)  followed  by  multiple  fiber-splitting  in  the  252~layer.  In  particu¬ 
lar,  there  is  no  fiber  breakage  involved  in  the  latter  case,  see  Figure  22. 

An  energy  analysis  of  the  respective  failure  modes  indicated  that  the 

transverse  cracking  mode  in  the  [+2S/+25]  is  suppressed  by  the  actual 

“  s 

thickness  of  the  25°-layer;  there  is  not  enough  strain  energy  to  cause  this 
failure  mode  at  a  lower  load  level.  Hence,  it  forces  a  fiber  breakage 
Instead.  However,  by  doubling  the  layer  thickness,  in  the  case  of  [252/— 25^ ] 
the  transverse  cracking  mode  emerges  at  a  much  lower  load.  The  difference 
in  the  failure  modes  results  in  the  difference  in  the  applied  ultimate  loads. 

Clearly,  the  strength  of  the  so-called  "load-carrying*'  plies  in  a  lami¬ 
nate  depends  also  on  the  particular  fracture  modes  that  cause  failure.  In 
the  examples  above,  the  ply-thickness  factor  is  seen  again  to  play  a 
profound  role  in  the  final  failure  mechanisms. 

When  the  respective  load-carrying  plies  are  laminated  with  a  90°-layer, 
such  as  in  [+25/+25/90-]  and  [25-/-25-/90™ ]  ,  transverse  cracks  in  the 

2  S  2  22S 

90°-layer  and  edge  delamination  are  found  to  proceed  any  failure  in  the 
25°-plies.  From  the  results  displayed  in  Table  A,  it  is  seen  that  the 
sub-laminate  cracking  events  neither  changed  the  failure  mode  of  the 
+25°-plies,  nor  decreased  the  final  ultimate  loading.  Furthermore,  the 
presence  of  the  90°-layer  in  [257/-257/90- ]  actually  increased  the  lami- 


nate  load-carrying  capacity.  The  reason  is  that  [252/-252]g  has  four  (4) 

25°-plies  in  the  core,  while  [252/-252/902]g  has  only  two  (2)  25°-plies 

scattered  by  the  presence  of  the  90°-layer  in  the  core. 

However,  adverse  effects  of  sub-laminate  cracking  on  laminate  final 

failure  modes  are  seen  from  the  test  results  displayed  in  Table  B.  Here, 

the  failure  modes  in  [+9/90-]  and  [+9/90o]  are  compared.  Note  that  the 

2  s  —  os 

90°-layer  thickness  in  the  latter  is  exaggerated  in  order  to  separate 

widely  the  respective  sub-laminate  crack  modes. 

In  the  cases  of  [0-/90.]  and  [+25/90n]  ,  the  final  failure  mode  in 
2  8s  —  8  s 

the  load-carrying  plies  involves  multiple  fiber-splitting,  which  occur 
along  the  root  of  a  90°-layer  transverse  crack;  see  Figure  22.  On  the 
other  hand,  there  is  no  fiber-splitting  in  the  [0./90  ]  or  the  [+25/90,] 

Z.  Z  3  Z  S 

laminates. 

A  three  dimensional  stress  analysis  will  show  that  there  is  a  stress 

concentration  at  the  90°-layer  transverse  crack  root;  both  the  applied 

tension  and  thermal  cooling  contribute  to  the  stress  concentration.  In 

particular,  the  size  of  the  stress  concentration  zone  is  proportional  to 

the  908-layer  thickness;  which  explains  why  the  considerable  difference 

in  the  final  loads  between  [0„/90_]  and  [0„/90_]  ;  and  between  [+25/90.] 

22s  2  8  S  —  2s 

and  [+25/90o]  ,  see  Table  B. 

—  8  s 

Fiber-splitting  seems  to  be  the  striggering  failure  mode  in  the 
[0./90o]  and  [+25/ 90_ ]  laminates.  The  splitting  mechanisms  is  due  to 

2  8  S  —  8  S 

the  existance  of  the  tensile  stress  normal  to  the  fibers,  which  is  magni¬ 
fied  in  amplitude  by  the  908-layer  transverse  crack  root. 

Similar  stress  condition  also  exists  in  the  [+30/90- ]  laminates.  But 

8  s 

in  this  case,  the  tensile  stress  (normal  to  +30°-f ibers)  is  not  large  enough 


to  cause  fiber  splitting;  see  also  Table  B.  As  a  result,  the  final  fialure 

modes  for  [+30/9CL]  and  [+30/90-]  all  involve  fiber  breaking;  and  their 
—  /  s  —  os 

final  failure  load  are  essentially  the  same.  Thus,  a  mere  increase  in  the 
angle  from  25°  to  30°  changes  the  failure  mode. 

Clearly,  the  final  failure  modes  in  the  "load-carrying"  plies  can  be 
influenced  by  nature  from  the  sub-laminate  crackings.  Moreover,  there  are 
probably  a  host  of  other  possible  modes  that  are  competing  for  dominance. 
The  specific  failure  occurrence  in  a  given  laminate  is  generally  determined 
by  the  lamination  geometry  as  well  as  the  intrinsic  material  properties. 
And,  of  course,  there  are  always  the  statistical  uncertainties  in  the  "in¬ 
trinsic"  material  properties. 

CONCLUSIONS 

This  paper  has  attempted  to  give  an  overview  of  a  fracture  mechanics 
approach  to  some  of  the  frequently  observed  sub-laminate  cracks  in  epoxy- 
based  composite  laminates.  The  laminates  discussed  in  the  paper  served 
mainly  to  illustrate  how  a  certain  simple  failure  sequence  is  developed 
and  how  the  predictive  model  is  constructed.  They  are  certainly  not  ones 
which  may  be  used  in  practical  design  applications.  The  main  objective 
is  to  enable  to  see  through  the  various  fracture  mechanisms  and  try  to 
understand  them;  at  times,  one  even  exaggerates  the  geometric  parameter 
to  a  radical  extent. 

Indeed,  the  field  of  fracture  in  composites  seems  infinitely  complex. 
Yet,  numerous  past  and  present  investigators  have  toiled  the  field  with 
many  successes.  And,  new  discoveries  are  being  continuously  made  to 
further  its  advancement.  The  results  reported  in  this  paper  represent 
but  one  effort  to  fill  the  dimensional  span  from  micromechanics  to  struc¬ 
tural  mechanics  in  composite  failure  analysis. 
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Fig.  12  Comparison  of  Experiment  and  the  Monte-Carlo  Simulations  of  Transverse  Cracks. 
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Fig.  20  SCHEMATIC  OF  THE  FRACTURE  SEQUENCE  IN  THE  (+25/90  )  LAMINATES 
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(a)  JUST  PRIOR  TO  EDGE  DELAMINATION,  (b)  SUBSEQUENT  TO  EDGE 
DELAMINATION,  (c)  JUST  PRIOR  fO  FINAL  FAILURE 
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Figure  22  Final  Failure  Modes:  (a)  Cross-sectional  Fiber 

Break;  (b)  25-layer  Transverse  Cracks;  (c)  0-layer 
Fiber  Splits  over  90-layer  Transverse  Cracks. 
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